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Abstract 

In the context of the seesaw mechanism, it is natural that the large solar and atmo- 
spheric neutrino mixing angles originate separately from large 2x2 mixings in the neutrino 
and charged-lepton sectors, respectively, and large mixing in the neutrino couplings is in 
turn more plausible if two of the heavy singlet neutrinos are nearly degenerate. We study 
the phenomenology of this scenario, calculating leptogenesis by solving numerically the set of 
coupled Boltzmann equations for out-of-equilibrium heavy singlet neutrino decays in the min- 
imal supersymmetric seesaw model. The near-degenerate neutrinos may weigh ^ 10 8 GeV, 
avoiding the cosmological gravitino problem. This scenario predicts that Br(fi — > ej) should 
be strongly suppressed, because of the small singlet neutrino masses, whilst Br(r — > /ry) 
may be large enough to be observable in B-factory or LHC experiments. If the light neu- 
trino masses are hierarchical, we predict that the neutrino less double-/? decay parameter 
m ee « jAm 2 sol sin 2 12 . 
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Neutrino oscillation data ||T], are converging towards unique solutions for both the solar 
and atmospheric neutrino anomalies |J. There are two large, almost maximal, mixing angles 
#12 and #23 in the light neutrino mass matrix, that give rise to the solar and atmospheric 
neutrino oscillations, respectively, whilst the third mixing angle # 13 is constrained to be 
small |5|. This pattern motivates theoretical approaches based on two 2x2 mixings 
instead of one general 3x3 mixing. 

The smallness of the neutrino masses is generally explained via the seesaw mechanism || , 
described by the superpotential []: 

W = N<{Y v ) i3 L 3 H 2 - Ei(Y e )ijLjHi + ^N c t (M N ) t ,N^ + iiE 2 E x . (1) 

Here the indices i,j run over three generations and M/v is the heavy singlet-neutrino mass 
matrix. We shall work in a basis where (M/v)ij is real and diagonal, (Mn)^ = M^Sij, and 
we define < M^ 2 < Mn 3 . In analogy with the CKM mixing matrix in the quark sector, 
the mixing matrix Vmns measured in neutrino oscillations is a product of two matrices 
Vmns — U\U U) where U u diagonalizes the light-neutrino seesaw mass matrix A4 U 

M v = Yj (Mn)- 1 » 2 sin 2 0, (2) 

according to 

U?M V U V = M D V , (3) 
and U e helps diagonalize the charged-lepton Yukawa coupling matrix Y e in ([!]): 

V^YJJe = Y e D . (4) 

Deriving the observed neutrino mixing angles from the neutrino Dirac Yukawa couplings 
(Y u )ij and Majorana masses Mn. is problematic. A generic difficulty, in view of the hierarchy 
^ m L/ ^ Am^ im , is that large mixing angles for both solar and atmospheric mixings can 
be obtained only at the price of some fine tuning. Technically speaking, the (23) sub- 
determinant of the light-neutrino mass matrix must vanish. While Yukawa textures which 
may accommodate this feature can be obtained from well-motivated physics ideas such 
as GUTs, Abelian and non-Abelian flavour symmetries, democratic principles, etc., many of 
these approaches feature unknown model coefficients of order unity that cannot be predicted 
without further assumptions, failing which they must be tuned a posteriori. 

1 We assume low-energy supersymmetry, which leaves unaltered the flavour parameters, while providing 
extra low-energy observables. 
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One way to attack this problem is to note that the observed leptonic mixing pattern can 
be regarded as the combination of two 2x2 mixings in the neutrino and the charged- lept on 
Yukawa couplings (Y v )ij and (Y e )ij, respectively. If the matrices U v and U e contain non-trivial 
mixings only in the (12) and (23) sub-matrices, respectively, large mixing angles #12, #23 and 
vanishing 6^ follow naturally. Notice, however, that the matrices U v and U e are of different 
nature. Whilst U e directly rotates the superpotential couplings in ([!]), U u diagonalizes the 
effective mass matrix obtained from the superpotential couplings via the seesaw relation (fj). 
Therefore, one may expect some differences between these two matrices U u and U e , which 
may be the reason why the solar mixing somewhat differs from the atmospheric one 0. 

Such a structure for Y v and U u would arise naturally if there is an (approximate) 5*2 
symmetry between the first and the second generation fields: Ni <-» N 2 and Li <-> L 2 , 
implying 

in the (12) sector. The large mixing angle 612 diagonalizing Y v follows immediately from 
this approximate symmetry, as does the prediction that Mjvj ~ M^ 2 , , while M^ 3 may 
be different. Larger flavour symmetries would be needed to extend this approach to the 
third generation, to incorporate the Ef fields, and to explain the largeness of #23- We do not 
pursue such model-building issues here, but rather pursue the phenomenological implications 
of the proposed structure (f|) for Y u , Y e and Mjv, namely that the Dirac Yukawa matrices 
Y v and Y e have non-trivial 2x2 (12) and (23) sub-matrices and non-vanishing diagonal 
elements (33) and (11), respectively, whilst the rest of their entries vanish and Mn has a 
pair of (almost) degenerate eigenvalues. We incorporate the low-energy neutrino data into 
our parametrization of the Yukawa couplings, and systematically scan over all the remaining 
free parameters of the seesaw model. In this way we include theoretical models of |7j that 
predict similar patterns. 

We focus our attention on leptogenesis |J in this framework, finding that the pair of near- 
degenerate heavy singlet neutrinos enable one to lower the scale of thermal leptogenesis below 
the gravitino bound on the reheating temperature of the Universe in supergravity and 
gauge-mediated models [|I(J. We study implications of this scenario for /3/3o„ decay |TiJ and 
on charged- lepton flavour- violating (LFV) decays . A important aspect of this scenario 
is that there are less free parameters than in the general seesaw model, implying testable 
phenomenological consequences. We find that r — > /ry may be observable, whilst /x — > is 

2 The above breakdown is basis-dependent: one could choose to work in the basis in which is diagonal, 
in which case the mixing in U e simply moves into Y v , while the physics remains unchanged. 
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suppressed, and make a specific prediction for neutrinoless double-/? decay. 

We start by counting the physical parameters of the model. In the chosen basis the 
Yukawa matrix Y v contains five complex parameters, and the basis for Ni is completely 
fixed. Since three phases can be removed by redefining the Lj fields, only two phases are 
physical. One can parametrize Y v as 

(Y v )ij = Z* k Y^X kj , (6) 

where the unitary matrices X and Z contain only (12) mixing, and hence only one mixing 
angle. The matrix X is real, while the matrix Z may be written as Z = P1ZP21 where Z 
is a real matrix. In this 'high-energy' basis, the diagonal matrices P\$ = diag(e 1,2 ,1, 1) 
contain the two physical CP- violating phases. Diagonalizing Y e according to (|4]), the real 
parameters correspond to three diagonal Yukawa couplings Y e D and one mixing angle in 
each of the rotation matrices U e and V e . Again, three out of five phases can be absorbed 
into redefinition of right-handed fields Ef. Thus the basis for E\ is now completely fixed 
and there is one physical phase in the mixing matrix U e = U e Ps, where U e is real and 
P3 = diag(l, 1, e 3 ). The mixing and one phase in V e are unobservable. This implies that, 
in the basis in which the charged lepton masses are diagonal, we have a total of 9 physical 
parameters in the neutrino Yukawa couplings, which together with the 3 unknown heavy 
masses make a total of 12 parameters in this version of the minimal seesaw model. This 
should be compared with the 18 physical parameters in the general case, implying more 
predictivity and hence better possibilities to test the scenario at high and low energies. 

We recall there are three types of leptonic observables in supersymmetric seesaw models: 
(i) leptogenesis via out-of-equilibrium decays of heavy singlet neutrinos and (ii) renormalization- 
induced mixings in the slepton mass matrix due to the off-diagonal Y v couplings, as well as 
(iii) the light neutrino masses and mixings. The leptogenesis CP asymmetry produced in 



the out-of-equilibrium decays of each of the Ni is given by EjL 14, 15 
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v^[Log(l + 1/xi) + , ( 7 ) 

where x\ = (M N J M Ni ) 2 . It is clear from (0) that the generated asymmetry depends only on 

Y v Yl = P*T(Y U D ) 2 Z T P 1 (8) 

and on the heavy neutrino masses. Notice that Y u Yj is independent of the left-rotations of Li, 
and is therefore independent of the Y e basis. Therefore, in our scenario the baryon asymmetry 
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of the Universe depends only on the single high-energy CP phase in Pi. If the CP-conserving 
parameters are known, this can be calculated from the observed baryon asymmetry. 

Renormalization of soft supersymmetry-breaking parameters due to the presence of Y v 
above the heavy-neutrino decoupling scales modifies the left-slepton mass matrix and 
trilinear soft supersymmetry-breaking A e terms. In the basis of diagonal charged leptons 



one has in leading-logarithmic order |L2 



(5m 2 
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which are proportional to 

y ] ly = p:uJxy d p 9 z t lz*p*y d x^Tj p p 



3- 



(9) 



(10) 



where L is a diagonal matrix: Ljj = \n(M GUT /M N .)Sij. Note that the high-energy CP- 
violating phases in (|T0|) are those in Pi and P3. The CP- violating observables in LFV 
processes all depend on a single CP- violating invariant J u = \mH12H23H31 HH, where 
H = YjLY v . This influences slepton physics at colliders and also determines the T-odd 
asymmetry in fi — > 3e P7| . Therefore, all low-energy CP- violating observables measure one 
combination of phases in P 2 3- Due to the near-degeneracy of the heavy neutrino masses, the 



lepton EDMs are unobservably small in this scenario [|T^ 

So far, we have considered the high-energy parametrization of the neutrino sector in 
terms of the diagonal Yukawa couplings and the unitary matrixes X, Z, U e . Counting 
physical degrees of freedom is straightforward in this approach, but it has limited appli- 
cability to low-energy phenomenology. Therefore we develop a complementary 'low-energy' 
parametrization. We start by discussing the light neutrino mass matrix @ in our scheme. 
Of course, at low energy it is most convenient to work in the basis in which charged leptons 
are diagonal. After redefinitions of the Li fields, the diagonalizing matrix for light neutrinos 
becomes 



V„P , 



and P = diag(l, e*^ 2 , e*^ 3 ) contains 



where V v is real and contains the mixing angles #12, #23 
two Majorana phases 2 ,3, whilst the third angle #13 vanishes and the neutrino oscillation 
phase 5 is absent. Therefore, all the low-energy neutrino observables, such as neutrino 
oscillations, /3/3o u decay, etc., depend on the 7 effective low-energy parameters, which are 
functions of the 12 parameters in ([!]). We recall that, whilst neutrino oscillations measure 
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the mass-squared differences of neutrinos and their mixing angles, (5(5q v decay measures one 
particular combination of their masses and mixing matrix elements: 



\vn, Pf ,\ = 



m u cos 2 9i2 + m U2 sin 2 9 12 e 



2i<t> 2 



(12) 



The effect of the Majorana phase 02 becomes visible only if m ui ~ Tn U2 , but its measurement 
is not straightforward even in this case |Tl|. In the case of hierarchical neutrino masses 



<C m U2 , we have the definite prediction m e 



sin 2 0i 



In order to satisfy automatically the oscillation data, the light neutrino parameters must 



be an input of the parametrization. We therefore define [19 



y/M^RJM^Ul 



(13) 



v sin (3 

where R is a complex orthogonal matrix. In our case, it has a non-trivial 2x2 (12) submatrix 
only, and is parametrized by just one complex number. Examining the combination Y u Yj 
using (0), we observe that the single complex phase in R is responsible for leptogenesis. On 
the other hand, @ tells that in this parametrization the renormalization-induced low-energy 
CP violation depends also on the Majorana phases in a non-trivial way. One can exploit 
this to take a complete bottom-up approach, since in the supersymmetric seesaw model the 
low-energy degrees of freedom may in principle be used to reconstruct all the high-energy 
neutrino parameters ]20| . A parametrization related to the solutions of the renormalization- 
group equations (RGEs) for the soft supersymmetry-breaking slepton masses (^) was worked 
out in g. 

Our central objective in this Letter is to study the phenomenology of leptons in the 
proposed scheme. In particular, we focus our attention on leptogenesis and its relations to 
neutrino masses, f3(3o u decay and LFV observables. A strong motivation is provided by the 
gravitino problem in generic supergravity theories, which restricts the maximum reheating 
temperature Tr of the Universe after inflation. In supergravity models with a gravitino mass 
of order 500 GeV or in some gauge-mediated models [|1(J the bound is 



T R < 10 8 GeV. 



(14) 



On the other hand, if the baryon asymmetry of the Universe is due to the decays of the 
lightest singlet neutrino N ly there is an M Nl - dependent upper bound on the CP asymmetry 
In the context of thermal leptogenesis, the observed baryon asymmetry implies a 



24 



ei 

lower bound <; 10 10 GeV J|5| which is potentially in serious conflict with (|i~4]). 



implications of particular neutrino mass textures on leptogenesis have been studied in Il4l 151 B2l [23| 
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This problem could be overcame by abandoning thermal leptogenesis and considering 
non-thermally produced neutrinos |26] , but such scenarios are still speculative and lack pre- 
dictivity. Another well-known solution to the problem in the context of thermal leptogenesis 
is to consider the decays of two heavy neutrinos which are approximately degenerate in 
IBf. Because the self-energy contribution to ej is enhanced for degenerate heavy neu- 



mass 



trino masses, the observed baryon asymmetry can be generated by moderately degenerate 
heavy neutrinos that are relatively light, which is the option pursued here. Since, in the 
proposed scenario, the Y v couplings spanning the (12) submatrix are separate and distinct 
from the (33) element, it is natural that the masses of the heavy neutrinos iVi^ differ from 
that of N 3 4 . 

In the minimal supersymmetric seesaw model, the baryon asymmetry originates from 
out-of-equilibrium decays of neutrinos and sneutrinos. Assuming M Na > M N2 ss M Nl , we 
need consider only the two lighter neutrino decays. Even if M^ 2 « , the CP asymmetries 
ei t 2 may be very different in magnitude. Although the products of the numerator and mass 
factors in (0) are identical for ei >2 , the denominators (Y u Yj) n and (Y u Yj) 2 2 m &y be very 
different f\. Therefore, Nip may contribute to the lepton asymmetry and to the reaction rates 
of the processes involved in very different ways, making the washout processes non-trivial. 
The dominant processes which determine the order of magnitude of the asymmetry are the 
AL = 1 neutrino and sneutrino decays. The AL = 2 scatterings are suppressed by additional 
powers of small Yukawa couplings, and are completely negligible in the low-M^ regime we 
consider here W^- There are also AL = 1 two-to-two scatterings involving top quarks 
and squarks. Those have additional powers of yf/ (Air) and make a contribution of relative 
order unity to the washout processes: we neglect them for simplicity. In supersymmetric 
models there are also processes transforming leptons into scalar leptons and vice versa, e.g., 
e + e <-> e + e , which we do take into account. 

In this approximation, defining Yj± = Y~ ± Y~\ for the scalar neutrinos and their an- 

3 j 

tiparticles, the Boltzmann equations for the ratios of particle densities divided by the entropy 
density, and for the lepton asymmetries Y L and Y Lg in fermions and scalars, respectively, 
are given by [ |15|| : 

dY N . , , 

3 - 1 3 - 1 7 



dz 




4 We assume here that _/Vi j2 are not exactly degenerate in mass, presumably because the (unspecified) 
underlying symmetry principle for first two generations, such as the S2 example mentioned earlier, is weakly 
broken. 

5 This would require the two-generation symmetry to be broken in the Yukawa couplings. 
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(15) 



where the temperature dependence is via z = M^/T, Y x eq denotes the equilibrium density of 
the particle denoted by the subscript x. The reaction densities for neutrinos and sneutrinos 
are given by [O 



In, = 2 % 



a iV a j 



Air 3 ^ " ' 33 z 
where a,j = (M/y^./M/yJ 2 , and that for the MSSM processes is 



KliZ-y/Oj 
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(16) 



(17) 



where m~ is the photino mass which contributes toe + e^e + ein the t channel. The total 
lepton asymmetry Yl = Y L + Y Lg is then converted by sphalerons into a baryon asymmetry: 
Y B = CY L} where C = —8/15 in the MSSM. 

In our subsequent numerical analysis, we fix the known light neutrino parameters to be 
Am 2 2 = 3 x 1CT 3 eV 2 , Am^ = 5 x 1CT 5 eV 2 , tan 2 9 2 3 = 1 and tan 2 12 = 0.4, corresponding 
to the LMA solution for the solar neutrino anomaly. The angle #13 is predicted to vanish, 
together with the neutrino oscillation phase 5, and the connection between leptogenesis and 
the oscillation phase proposed in [23[] does not hold. We assume the normal mass ordering 
for light neutrinos, m Vl < m U2 < m U3 , and generate Y v according to (0). All the unknown 
input parameters are generated randomly as follows. The lightest light neutrino mass m Ul 
is generated in the range (10~ 5 — 1) eV, the Majorana phases 02,3 in the range (0 — 2ir), the 
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lightest heavy neutrino mass in the range (10 6 - just above 10 8 ) GeV, the degeneracy 
parameter A = ln(Mjv 2 / — 1) in the range (—15 — 1), the heaviest neutrino mass 
in the range (2 x Mn 2 — 10 15 ) GeV and the complex parameter tyi of the orthogonal matrix 
R in the range 1 7" 12 1 = (0 — 10) with an arbitrary phase. Each mass parameter is generated 
with a flat distribution on a logarithmic scale. We solve the Boltzmann equations ( |I5D 
numerically, assuming initial thermal abundances for neutrinos and sneutrinos, and require 
that the induced baryon asymmetry be in the range 2g 



3 x 10" 11 & Y B <, 9 x 10" 11 . (18) 

We require that the neutrino decay width satisfy < 10 -5 (Mjv 2 — Mzvi), so that we are far 
from the resonant region and can trust our perturbative calculation. Subsequently we study 
correlations between the leptogenesis parameters and the light and heavy neutrino masses, 
the PP„o decay parameter m ee , and the LFV decays of the charged leptons. 



o 




7 5 -5 -2.5 

Ln(M N7 /M N1 -l) 



Figure 1: Scatter plot of the lightest singlet neutrino mass as a function of the degeneracy 
parameter • In (M^/M^ — 1). The baryon asymmetry is required to be in the range (fl6[). 

We find a weak correlation of the allowed range of M Nl with A = \n(M N2 /M Nl — l), as seen 
in Fig. [I]. We see that successful leptogenesis is possible with light A^ 12 already if they are 
degenerate in mass at the level of a few percent. Motivated by the gravitino problem, we limit 
MjVi to the range 10 6 GeV to about 10 s GeV in our subsequent numerical examples. We find 
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a stronger correlation between the leptogenesis parameters and the neutrino degeneracy 
parameter A. In Fig. |2] (a) we plot the CP asymmetries ei^, denoted by black and green 
points, respectively, as functions of the degeneracy parameter A = ln(Mjv 8 /MjVj — 1). There 
is a A-dependent upper bound on ei )2 , coming from the mass function in (0). Solutions to 
the Boltzmann equations (fT5|) are often discussed in the literature in terms of effective mass 
parameters 



ITLi = 



(w?)..^. da) 



If leptogenesis originates from N\ decays, as is the case for hierarchical heavy neutrinos, 
fixing Yb implies almost a one-to-one correspondence between e\ and mi in the low 
region |T3[ |25f . However, the solutions to (|T5|) are non-trivial, as seen in Fig. (b), where 
we plot £x t 2 versus the parameters respectively. A direct correlation is observed only 
for high values of fh\ ~ m 2 ^ 5 x 10 -2 eV, where both neutrinos AT 1)2 contribute to Yb in 
a similar way. This region corresponds to large and quasi-degenerate light neutrino masses. 
Because of strong washout effects, ei j2 must be large and, correspondingly, M^/M^ — 1 
small in this region. For lower values of rh\ the non-linear shape of Fig. |2| (b) indicates the 
non-trivial washout effects of A^ 2 . 

To study implications of our neutrino mixing scenario and leptogenesis for low-energy 
observables, we plot in Fig. |] (a) the lightest neutrino mass m ui versus the (3(3q v parameter 
m ee . There is an upper bound on m Vl coming from the requirement of successful leptogenesis. 
This is an artifact of our lower bound on the degeneracy parameter A and is not physical. 
If we allow smaller A, £i )2 are enhanced and the neutrino mass m Vl can be higher. At 
high values of m Ul , there can be cancellations in m ee flT2|) , because the phases </> 2 ,3 are free 
parameters not correlated with leptogenesis. However, for small m Vl there is the definite 



prediction m ee = yAm^j sin 2 12 , which is below the sensitivity of the currently proposed 
j3j3 0u decay experiments |29|] Q. 

To study LFV in our scenario, we fix the soft supersymmetry-breaking parameters at the 
GUT scale to coincide with one of the post-LEP benchmark points |30]]: mi/2 = 300 GeV, 
m = 100 GeV, A = 0, tan/3 = 10 and sign(fi) = +1. We plot in Fig. | (b) the 
renormalization-induced branching ratio of the decay ji — > ej versus the branching ratio 
of r — >• /ry. Because we require M Nl to be relatively low, the Yukawa couplings Y u related to 
first two families are also very small, suppressing Br(fi — ► ej) below the presently observable 
level. Since the branching ratios scale as tan 2 /3, reaching the level 10~ 14 or 10 -15 planned to 



6 If 9i 3 is not exactly vanishing, there is an additional contribution to (12) and the prediction for m e 
becomes less certain fUJI. 
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Figure 2: Scatter plots of the CP-violating asymmetries ei >2 denoted by black and green (grey), 
respectively, as functions of the degeneracy parameter \n(M^ 2 / M^ — l) and the effective mass 
parameters The baryon asymmetry is required to be in the range ftTBi). 




Figure 3: Scatter plots of the lightest neutrino mass m Vl versus the (3j3Q U -decay parameter 
m ee , and Br(fM — > ej) versus Br{r — > fij). The baryon asymmetry is again required to be in 
the range (fJ^j. 
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be achieved at PSI [3I| and the neutrino factory experiments P2"| |, respectively, is unlikely for 
Mn 1 ~ Mtv 2 ^ 10 s GeV. This indicates that the solution to the gravitino problem and the 
non-observation of 11 — > in present experiments may be related. However, Br[r — *> fij) 
is not necessarily suppressed, because N 3 can be heavy and the related Yukawa couplings 
large. Therefore r — > /ry may be observable in B-factory or LHC experiments, which should 
achieve a sensitivity Br{r — > fj,j) ~ to 10~ 9 [33| . 

In conclusion: we have studied a scenario in which the large leptonic mixing angles #12 
and 823 originate from 2x2 mixings in the neutrino and the charged lepton Yukawa matrices, 
respectively. Using a convenient phenomenological parametrization and scanning over the 
all free seesaw parameters, and motivated by the gravitino problem, we have focussed on 
lowering the thermal leptogenesis scale by postulating moderate degeneracy of two lightest 
heavy neutrinos. This scenario turns out to be quite predictive for j3j3 0u decay and for the 
LFV decays of charged leptons, since it involves less free parameters than the general seesaw 



model. We predict suppressed fx — > &y, m ee = \J Am sol sin 9\2 for hierarchical neutrinos and, 
in general, an observable rate for r — > ii r ). 
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